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Rocket observations at high latitudes indicate a 
significantly cold thermal structure during summer in 
the upper mesosphere and lower thermosphere. From de- 
tailed considerations of radiative transfer we have ob- 
tained the infrared heating rates in this region of the 
atmosphere and deduced vertical velocities from these 
heating rates. 

We find that strong infrared flux convergence in the 
vicinity of 85 km results in a heating of about 10° K/day 
and positive vertical velocities of about 5 cm/sec. Such 
vertical velocities favor the formation of noctilucent 
clouds at this altitude. 
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Vertical Motions in the Upper Atmosphere Deduced 
from Rocket Measured Temperatures 


Introduction 

Recent high-altitude rocket-grenade experiments at Church- 
hill, Canada, Barrow, Alaska and Kronogard, Sweden have re- 
vealed some interesting features of the temperature structure 
in the earth’s mesophere and lower thermosphere. These ob- 
servations indicate that, at high latitudes, this region of 
the atmosphere is significantly warmer in winter than in summer, 
as shown in Figure I (Theon, Nordberg and Smith, 1966) . Such 
a temperature pattern cannot be explained by considering a 
radiation balance alone. It has been suggested that, in 
winter, dynamical heating (Leovy, 1964) , chemical heating 
(Kellogg, 1961) and heating by dissipating gravity waves 
(Hines, 1963; Leovy, 1966) may play a role in establishing the 
thermal structure at these levels. In summer, however, dy- 
namical cooling appears to be the only significant non-radiative 
process involved. 

An examination of radiative heating and cooling in the 
stratosphere and mesosphere has been made by Murgatroyd and 
Goody (1958) using the information regarding the thermal struc- 
ture of the atmosphere available at that time. On the basis 
of their heating rates, Murgatroyd and Singleton (1961) de- 
rived a meridional circulation in this region of the atmosphere. 
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The extremely cold temperature structure in the vicinity 
of the mesopause in summer, which has been revealed by the 
recent observations prompts us to re-examine the radiative 
balance of this region. 

For the observed summer temperature structure at Kronogard 
and Barrow (Figures II and III) the infrared heating of the at- 
mosphere at these levels is obtained from a detailed radiative 
transfer calculation in the 15 p C0 2 and 9.6 p 0^ bands. From 
this heating we derive the vertical velocities necessary to 
maintain the observed temperature structure, assuming solar 
ultraviolet heating is absent. Such estimates of vertical 
velocities are conservative , in our view, since any solar heating 
would only enhance vertical motions. 
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Theoretical Considerations 


From the theory of radiative transfer, it is possible to 
calculate an equilibrium temperature structure when the sources 
and sinks of radiative energy within the atmosphere are known. 
Conversely, if a thermal structure is specified, for example, 
from observations, and if an assumption of quasi-equilibrium 
is made, the radiative heating and cooling rates in the at- 
mosphere can be obtained. 

The basic equation of radiative transfer relates the net 
flux divergence to the absorbed and emitted energies per unit 
volume and per unit time: 


dl (V, 8 , z) <jv duu = - JJ I (v, 0 , z) n (z)k (v, z) dv duu 
ds 

+ J(v, z) n (z)k (v, z) dv duu 

where the symbols have the following meaning: 
v frequency of radiation 

w solid angle 

9 zenith angle 

s geometric path length (ds = sec 9 dz) 

z altitude 

n(z) number density of absorbing molecules 

k(v,z) absorption coefficient 

J(v,z) source function 


(1) 


I(v,9,z) intensity of radiation , 


- 3 - 


With the boundary conditions of zero downward flux of 


infrared energy at the top of the atmosphere z„, and a black- 
body surface at the ground z Q ,the basic equation of radiative 
transfer (1) can be solved to yield at any level z^, the up- 
ward and downward intensities, lf(v,0, z^) and 
respectively: 

Z 1 

St (v , 0 , z 1 - z) 

w 5z 
z o 

Z H 

St ( v, 0 , z - -1 )^ 

Sz 


= B (v, z q ) t (v, 0 , - z q ) + J J ( v , z ) - 


H 

ll(v,0, Zl ) = J J(V,z) 


dz (2) 


(3) 


B(v,z q ) is the Planck intensity from the ground, and t(v,0,|z, - z|) 
represents the transmission of the atmosphere at frequency v 
and zenith angle 0 between levels z and z^: 


t(v,0,|z 1 - z | ) = 1 - A(v,0,|z 1 - z | ) (4) 

where A(v,0,|z^ - z|) is the fractional absorption. 

From expressions (1), (2) and (3) the intensities and thereby 

the net flux divergence can be evaluated, provided the trans- 
mission function t and the source function J are known. A dis- 
cussion of the transmission function is given in the next section. 

The source function J is commonly assumed to be the Plank 
intensity B. However, such an assumption fails in the 15 p CC^ 
band at low pressures in the high atmosphere, where molecular 
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collisions are not frequent enough to maintain the vibrational 
energy levels in a Boltzmann distribution. According to Curtis 
and Goody (1956) the source function in this band when there is 
departure from local thermodynamic equilibrium is given by: 

J(v 

J J 

where cp is the radiative lifetime and X is the vibrational re- 
laxation lifetime. 

Since the source function J shown in equation (5) is itself 
dependent on the flux divergence, ^ ^ds^ ~ ~ ~ in the 15 p 

band, an iterative procedure is adopted to determine J. 

Radiative transfer in the 9.6 p 0 3 band should be treated 
in a similar manner to allow for possible departure from thermo- 
dynamic equilibrium. A lack of knowledge of the relaxation 
time for the 0^ molecule prevents us from doing so, and we have 
assumed the validity of Kirchhoff 's law in the 9.6 p band. At 
high altitudes the flux divergence of energy in the 9.6 p band 
is small compared to that in 15 p CC> 2 band, so that such an 
assumption does not lead to significant error. 

Taking the temperature structure in the atmosphere from 
rocket observations we have computed the total infrared flux 
divergence in the 15 p C0 2 band and 9.6 p 0 3 bands. Under steady 
state conditions, this total flux divergence must be balanced 
by solar energy absorption and dynamical heating, so that 



- 5 - 


(6) 


/ 9T(z) g 
-W(z)p(*)C ( - 

p 

where W(z) is vertical velocity 
p (z) is mass density 

C is specific heat at constant pressure 
P 

g is the acceleration due to gravity 
Q(z) is the solar energy absorbed per unit volume and 
per unit time by 0^ and 0^ . 

Equation (6) can be solved for W(z) provided the other 
quantities involved are specified. 

The solar ultraviolet energy absorbed by 0 2 and 0^ can 
be evaluated by specifying the absorption cross sections, spectral 
energy distribution, and the local concentrations of 0 2 and 0^. 

The vertical distribution of these gases may be obtained from 
photochemical considerations (See, for example, Prabhakara and 
Hogan, 1965) . 


) + Q(Z) = jj ' n p -L l'*) dvdw 
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Transmission Functions 


The relative importance of Lorentz and Doppler broadening 
depends upon the region of the atmosphere under consideration. 

In the optically thick layers of the lower atmosphere, where 
collisions are frequent, line broadening is predominantly Lorentz 
in character. At very high levels, on the other hand, the Doppler 
broadening is important, while at intermediate heights, a combin- 
ation of both Lorentz and Doppler effects is present. A general 
expression for the variation of the absorption coefficient over 
a spectral line which takes both phenomena into account is given 
by 

k(V) = S j ( /ttcci d ) exp " ( c a D ) tt (v - V q + mTTcp + (7) 

Here, the symbols have the meaning: 

S total line intensity 

U molecular velocity 

v q central frequency of line 

c speed of light 

a Lorentz half-width of line 

J-i 

a D Doppler half-width of line. 

For small values of the parameter a = a^/a^, Harris (1948) 
has expressed this integral as a power series in a, with co- 
efficients H^ which are tabulated as functions of the parameter 
v = (v - v Q )/ a D ’ 
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( 8 ) 


k = —■ — C H o (v) + aH i (v) + a 2 H (v) + . . .] 
% / TT a D ol2 J 


This formulation was used to evaluate k v for a < 0.5. Using 
this expression, the mean fractional absorption for a spectral 
interval is 


A 



+d/2a 

I 

-d/2a D 


[1 


exp (-k v u) ]dv 


( 9 ) 


where d is the mean line spacing and u is the optical path. 
Equation (9) can be solved numerically by substituting the tab- 
ulated functions in equation (8) for k v - 

Van der Held (1931) has shown that for values of ct > 1, the 
mean fractional absorption for a spectral interval is closely 
approximated by the Ladenburg-Reiche (1913) formula: 


A = 1 - exp [-3xe _x (I Q (x) + I 1 (x))] (10) 

where 

2nO-L 

P = ~d 

Su 


2na L 


and I and I. are the Bessel functions of order zero and one 
o 1 

and of imaginary argument. 

The two equations (9) and (10) give the fractional absorp- 
tion due to one spectral line. When broad spectral intervals 
consisting of a large number of lines are considered, the trans- 
mission and absorption of a gas can be represented by any one 
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of several band models. Here we have adopted a "statistical" 
band model in which the mean line spacing is held constant 
throughout the band, while the line intensity is allowed to 
vary from interval to interval within the band. So, if the 
mean line spacing d for the whole band and mean line intensity 
S for each one of the spectral intervals within the band are 
known, the average transmission or absorption in each of the 
spectral intervals can be calculated from equations (9) and (10) . 

The 15 g CO 2 band was divided into six spectral intervals 
of 50 cm 1 width from 550 cm "*" to 850 cm ^ . The Lorentz half- 
width a L of the lines in this band is given as 0.064 cm”'*' by 
Kaplan and Eggers (1956) for a temperature of 298° K and a 
pressure of one atmosphere. Yamamoto and Sasamori (1958) have 
used a mean line spacing d of 1.55 cm ^ in their calculation of 
the absorption of this band. Thus, a value of the parameter p 
is determined as 0.24. With this knowledge of the magnitude 
of 0 , another parameter, S/2 tt 0. l , was deduced for the various 
spectral intervals from the transmission tables of Stull, Wyatt 
and Plass (1963). This parameter is given in Table IA. 

The 9.6 g 0^ band was treated as a single spectral interval 

extending from 990 cm ^ to 1090 cm ^ . Plass (1960) had made 

an analysis of the 9.6 g 0^ band area measurements by Walshaw 

(1957) and so had evaluated the parameters 0 and S/2ttoc t for this 

Xj 

band. It was found that the values obtained by Plass can re- 
produce the band area as measured by Walshaw quite accurately 
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and so they were adopted in our investigation. The Loren tz 
half-width of the spectral lines in the 9.6 u 0^ band is given 
by Goody (1964) as 0.076 cm ^ at 293° K and 1000 mb. Using this 
value, the quantities S and d can be obtained (Table IB) . 

The value of the Doppler half-width at a temperature 
of 225° K is also given in these tables for both bands. 

Using these values of the various band parameters, the 
fractional absorption is calculated from equations (9) and (10) 
for different path lengths and pressures. As an illustration, 
the absorption obtained for the 15 u CC^ and 9.6 u 0^ bands, 
plotted as a function of 2rrax, is shown in Figures IV and V. 

In these figures, lines of constant values of the parameter a 
are shown. Each of these lines corresponds to a line of con- 
stant pressure if the temperature remains unchanged at 225° K. 

For low pressures (small a) the absorption is mainly due to the 
Doppler broadening and, with the increase of path length, the 
absorption increases very slowly until the wings of the absorp- 
tion lines gradually begin to absorb. Then the absorption in- 
creases in square root fashion with path length until the level 
of saturation is reached. 

In calculation of atmospheric transmissions over long path 
lengths, the variation of pressure and temperature must be 
taken into account. Therefore, the weighting procedure de- 
veloped by Curtis (1952) and Godson (1953) was used to obtain 
the transmission. The temperature dependence of the line strengths 
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in the 15 u C0 2 band is treated in a way similar to that used by 
Praphakara and Hogan (1965) . The variation of a with temperature 
is also incorporated in the calculation. 
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Band Parameters 
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Results and Discussion 


We have examined the infrared radiative heating corresponding 
to the rocket-measured temperature profiles at Kronogard (66° N, 

29 July, 1963) and Barrow (71° N, 9 August, 1965) and determined 
the vertical velocities in these cases by the procedure outlined 
above. The heating rates and resultant vertical velocities are 
given in Figures II and III. 

The convergence of infrared flux between 70 km and 90 km 
altitude leads to a significant diabatic heating of the order 
of 10° K/day. This heating alone implies vertical velocities 
of ~ 5 cm/sec near 80 km. This is a consequence of a nearly 
adiabatic temperature gradient in this region of the atmosphere. 

Calculations of infrared heating rates were made for CC> 2 
vibrational relaxation times X of 10 ^ sec and 5 x 10 6 sec 
(at a pressure of 1 atmosphere) to illustrate the effect of this 
parameter on the vertical velocities. The vibrational relax- 
ation time of the fundamental of C0 2 (15 u band) is known to 
within an order of magnitude, 10 ^ sec to 10 ^ sec (Marriott, 

1964) . This relaxation time is inversely proportional to pres- 
sure and many have a weak temperature dependence (Herzfeld and 
Litovitz, 1959). We have considered only the pressure depen- 
dence as variations in pressure in the atmosphere are much larger 
than temperature variations. It can be seen in Figures II and 
III, that decreasing the relaxation time increases the infrared 
heating. Were the atmosphere in a state of local thermodynamic 
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equilibrium at these levels, heating rates in excess of 20° K/day 
and vertical velocities in excess of 10 cm/sec would be obtained. 

Infrared transfer of energy in the 6.3 u and far infrared 
rotation bands of H 2 0 was not taken into account in this study. 
Estimates of the water vapor content in the high atmosphere 
differ markedly (see Theon, Nordberg and Smith, 1966) . We be- 
lieve that the presence of significant amounts of water vapor 
would only augment the infrared flux convergence in the cold 
region of the atmosphere around 85 km, and result in larger 
vertical velocities. 

Solar ultraviolet heating by 0^ has not been considered 
here, as the interaction of 0^ with H 2 0 at these levels would 
reduce the concentration of 0^ and its heating effect quite dras- 
tically (Hunt, 1966; Hessvedt, 1965) . Since the amount of H 2 <D 
present near the mesopause is uncertain, elaborate photochemical 
calculations to determine the 0^ concentration were not made. 

Any ozone heating in the vicinity of the mesopause would contribute 
to increase the vertical velocities. 

In view of the uncertainties of the heating by H 2 0 in the 
infrared and by 0^ in the ultraviolet, a detailed determination 
of the meridional circulation involving horizontal heat transport 
has not been made here. The cold air advection which one might 
expect at higher levels in the summer mesosphere would tend to 
compensate in part for the strong infrared flux convergence at 
these levels. However, we feel that these estimates of vertical 
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velocities of ~ 5 cm/sec near 80 km, although somewhat larger 
than those of Murgatroyd and Singleton (1961) , ~ 1 cm/sec, are 
quite reasonable in view of the extremely cold temperatures ob- 
served at these levels. 

The vertical velocities derived here might be related to 
the occurrence of noctilucent clouds in the vicinity of the 
earth’s mesopause. An extensive review of noctilucent cloud studies 
is given by Fogle and Haurwitz (1966) . Noctilucent cloud dis- 
plays are observed to be quite widespread and thus imply the 
existence of large scale vertical motions. Charleson (1965) con- 
cludes that vertical velocities of the order of 20 cm/sec are re- 
quired to sustain a cloud of spherically shaped particles of 0.1 u 
size with a mean density of ~ 1.0 g/cc. On the other hand. Chap- 
man and Kendall (1965) estimate that, if the particles are of 
needle or flaky form with a mean density of ~ 0.1 g/cc, vertical 
velocities of ~ 1 cm/sec (as suggested by Murgatroyd and Singleton 
(1961)) would be sufficient for cloud formation. The vertical 
velocities of ~ 5 cm/sec which are suggested by this study, when 
extremely cold temperatures are present, would favor theories of 
noctilucent cloud formation based upon less restrictive assumptions 
with regard to particle shape, size and density. 
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FIGURE CAPTIONS 


Figure I : Typical high latitude summertime and wintertime 

temperature profiles. These averages include soundings from 
Churchill, Barrow and Kronogard. (After Theon, Nordberg and 
Smith, 19o6.) 

Figure II : Observed summer temperatures (heavy solid line) , 

calculated infrared heating rates and resultant vertical velocities 
((A) light solid lines, = 10"5 sec; (B) dashed lines, X c = 

5 x 10"° sec) at Kronogard, Sweden. 

Figure III : Observed summer temperatures (heavy solid line), 

calculated infrared heating rates and resultant vertical velocities 
((A) light solid lines, X 0 = 10~5 sec; (B) dashed lines, Xc = 

5 x 10"* sec) at Barrow, Alaska, 

Figure IV : Absorption by mixed Doppler and Lorents lines in 

the 15^ c0 2 band. The abscissa 2v<sx = Su/x^, , and d/d,,** 2.46x103. 

Figure V : Absorbtion by mixed Doppler and Lorentz lines in 

the 9.6/1 ^3 band . The abscissa 2-nax = Su/ct-e, , and d/ct** 66.3. 


o 

O') 


o 

00 


o o 

<x> 


(IAIX) 1H0I3H 


o o 

LO ^ 


M 


Figu 
















